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On the Early Time X-ray Spectra of Swift Afterglows I: 
Evidence for Anomalous Soft X-ray Emission 
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^ ■ ABSTRACT 

> 

O 

1^ . We have conducted a thorough and bhnd search for emission hues in > 70 

f^ . Swift X-ray afterglows of total exposure ~ lO^s. We find that most afterglows 

^ . are consistent with pure power-laws plus extinction. Significant outliers to the 

^ . population exist at the 5-10% level and have anomalously soft, possibly thermal 

OO . spectra. Four bursts are singled out via possible detections of 2-5 lines: GRBs 

^ : 060218, 060202, 050822, and 050714B. Alternatively, a blackbody model with 

O , kT ~ 0.1 — 0.5 keV can describe the soft emission in each afterglow. The most 

Q . significant soft component detections in the full data set of ~ 2000 spectra cor- 

respond to GRB 060218/SN 2006aj, with line significances ranging up to ~ 20a. 
A thermal plasma model fit to the data indicates that the flux is primarily due to 



Oh' 



P ! L-shell transitions of Fe at ~ solar abundance. We associate (> 4cr significant) 



H 



c/3 , line triggers in the 3 other events with K-shell transitions in light metals. We 

favor a model where the possible line emission in these afterglows arises from 



X 



the mildly relativistic cocoon of matter surrounding the GRB jet as it penetrates 
and exits the surface of the progenitor star. The emitting material in each burst 
p^ ', is at a similar distance ~ 10^^ — 10^^ cm, a similar density ~ 10^^ cm~'^, and 

subject to a similar flux of ionizing radiation. The lines may correlate with the 
X-ray flaring. For the blackbody interpretation, the soft flux may arise from 
break out of the GRB shock or plasma cocoon from the prog enitor stellar wind. 



as recently suggested for GRB 060218 (ICampana et al.ll2006l ). Due to the low z 



of GRB 060218, bursts faint in 7— rays with fluxes dominated by this soft X-ray 
component could outnumber classical GRBs 100-1. 

Subject headings: gamma rays: bursts — supernovae: general — X-rays: general 
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Introduction 



One of the key open questions in the study of Gamma-ray bursts (GRBs) is that of 
the X-ray afterglow hues. Claims of low to moderate significance emission lines have been 
made based on data from severa l missions: Fe lines have been detected in afterglow data 
from AS CA dVoshida et al.lll999h . Beppo-SAX dPiro et aljEoogl : kntoneUi et"aDl2002h . and 
Chandra ( Piro et al.l I2OOOI ): lines from highly ionized light, multiple-a elements like Mg , 



Si, S, Ar, and Ca have been dete cted in afterglow d ata from XMM ( JReeves et al.l 12002 



Watson et al.l 120031 ) and Chandra ( JButler et al.ll2003l ). The detections are challenging to 
explain because th ey typically imply large, concentrated masses of metals in the circumburst 
material (see, e.g.. lLazzati. Campana. fc Ghisellinilll999l ) and a very efficient reprocessing of 



the non-thermal afterglow continuum into line radiation (see, 
2OOII : iLazzati. Ramirez- Ruiz, fc Reesll2002l ). 



Ballantyne &: Ramirez-Ruiz 



Rutledge &: Sakd (120031 ) and ISako. Harrison, fc Rutledgd (120051 ) argue that the claims 



made to date lack the necessary significance needed to prove that the X-ray lines are real. 



We ad dress the line significance in the case of GRB 011211 in a recent paper (JButler. et al. 



2005al ) and find that the differing significance estimates are simply due to different input as- 
sumption in th e continuum modelin g. T o prove consistency betw een the analysis techniques 
utilized by the iReeves et al.l (120021 ) and iRutledge fc Sakd (|2003[ ) , we developed software to 
autonomously detect one or multiple emission lines. Below, we employ the line search tools 
to comb the vast Swift X-ray Telescope (XRT) data set. With nearly 0.4 years in total 
exposure for ^ 70 bursts observed typically within 1 hour of the GRB, the mounting XRT 
sample represents a unique laboratory in which we can test the veracity of historical line 
claims and probe the physics of new phenomena (e.g., flares, rapid and unusually flat light 
curves) at early times. 

As we discuss below, divided up in time, the majority (~ 90%) of Swift X-ray afterglows 
are well modeled by simple power-laws in energy, partially absorbed below ~ 1 keV by gas 
along the line of sight. The typical power-law photon index is F ~ 2. However, a small 
fraction (~ 10%) of the afterglows in the sample exhibit very different spectra (F ~ 5 power- 
laws or prominent residuals near 1 keV), which trigger the line detection robot. These 
spectra can be fit with a blackbody model in addition to a pow er- law, or with eraission lines 
in addition to the power-law. One such event — GRB 060218 (ICusumano et al.ll2006al) — is 



clearly associated with a SN 2006aj (JMirabal et al.ll2006l : ISoUerman et al. 



2006 



Modjaz et al. 



20061 ) and may be the smoking gun tying the anomalous soft X-ray spectra from this and 3 
other GRB afterglows to supernovae and possibly directly to the breakout of the GRB shock 
from the progenitor star. 
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Data Reduction and Continuum Fits 



Our automated pipeline at Berkeley downloads the XRT data in near real time. The 
burst right ascension and declination are gleaned from the Gamma-ray Coordinates Network 
(GCN) notices in order to run the xrtpipeline reduction task from the HEAsoft G.Qj 
software release. We use version 007 of the response matrices. From there, we bin the 
data in time, exclude pileup chip regions for each time interval, refine the X-ray position, 
and produce spectra using custom IDL scripts. For the XRT position counting (PC) mode 
pileup rejection, we discard regions with count rates in excess of 0.5 cts/s (0.5 - 10 keV) (e.g.. 



Nousek et al.ll2006l ). For the Windowed Timing (WT) mode data the cutoff is 150 cts/s. For 
each time slice, we determined whether or not this threshold has been surpassed by analyzing 
counts within a 3.5 pixel (inner radius) to 16 pixel (outer radius) annulus around the source. 
The inner radius is then expanded or contracted to bring the count rate close to (but not 
beyond) the pileup threshold. For ~ 99% of the data by time (or ~ 60% by mass), there is 
no pileup and the resulting inner radius is zero. 

Spectral response files are generated using the xrtmkarf task for each time slice, and 
the time slices are weighted by exposure and summed. The spectra are fit in ISIqj. For each 
spectral bin, we require a S/N of 3.5. We define S/N as the background-subtracted number 
of counts divided by the square root of the sum of the signal counts and the variance in the 
background. We define the background region as that where the number of counts in an 
aperture the size of the source extraction region is within 2-sigma of the median background 
over the chip in that aperture for one contiguous follow-up observation. In most cases, the 
background is entirely negligible [^ 1% of the source flux). 

We fit the integrated spectra for each afterglow and also slice the data in time to search 
for transient emission features. For the PC mode data, we consider data in the 0.3-10 keV 
energy range, while we restrict to 0.5-10 keV for the WT mode data, where the low energy 
response appears to be poorly calibrated (see Section [53] for more details). For the period 
between GRBs 050124 and 060313, we reduce data from 72 PC-mode and 44 WT-mode 
spectra, and sub-divide the PC and WT mode spectra into 208 and 944 spectra, respectively, 
each with ~ 500 counts. The time coverage of the time-sliced spectra are plotted in Figure 

m 



^http://heasarc. gsfc.nasa.gov/docs/software/lheasoft/ 
2http://space.mit.edu/CXC/ISIS/ 
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Fig. 1. — (A) Start (black) and stop (grey) times for the WT mode spectra used in our analysis. 
(B) Start and stop times for the PC mode data. The duration of the WT mode spectra range from 
3.5 s to 4.5 ksec, with a median duration of 17s. The PC mode spectra durations range from 113 
s to 2.4 Msec, with a median duration of 9.9 ksec. 



We test the XRT data against two models for the continuum emission: a power-law 
and a blackbody model, both with photoelectric absorption. Overwhelmingly, the XRT data 
are well fit by absorbed power-laws. There is a clear clustering of the photon indices F: 
the median and median absolute deviation about the median are F = 1.9 ± 0.3. (Figure 



|2B)- The median column density in excess of the Galactic value (JDickey &: LockmanlllQQOl ) 



is Nh — a?"// Galactic = (1.9 ±0.9) X 10^^ cm~^, marginally consistent with zero. The integrated 
(i.e., not time sliced) spectra show consistent values, F = 1.9 ± 0.2 and Nh — iV/^^caiactic = 
(0.9 ± 1.2) X 10^^ cm~^. Soft outliers to these trends are labeled in Figure [2] and discussed 
more below. The fluxes of the spectra cover a broad range (Figure EJG) . They are roughl y 
consistent with fluxes previously measured for afterglows by Beppo-SAX f lCosta et al.lll999l ). 
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Fig. 2. — (A) Reduced x^ values for the continuum fits, shifted to a fiducial number of degrees 
of freedom z^ = 34. The vertical line present the 90% confidence rejection threshold, to the right 
of which the fit should be rejected. A small fraction of the spectra (~ 10%) are acceptably fit 
with a blackbody model (grey line), however most (~ 90%) are well fit by absorbed power-laws 
(black line). (B) Continuum parameters Nh and photon index T for the absorbed power-law model 
fits. There are soft outliers (GRBs 050714B, 050822). (C) Flux vs. time interval plot for t he 500 
count spectra. Also, plotted is the mean afterglow flux from Beppo-SAX (jCosta et al.lll999l ). The 
Beppo-SAX X-ray fluxes are within ±1 dex of the red line. 



3. The Emission Line Search: Efficiency and False Detection Rate 

To search for lines in each spectrum, we fit unresolved emission lines in addition to the 
power-law continuu m in the 0.3-5.0 keV band. We allow for the possibility of multiple emis- 



sion lines (see, e.g., iReeves et al.ll2002l ). The lines are fitted successively, starting with one 
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line and adding a total of five lines. The best-fit line location at each ste p is found by con- 
volving the fit residuals with the Line Response Function (LRF) (see, e.g.. lRutledge fc Sako 
20031 ). The convolution tests for line energi es on the 5E = 5 eV A ncillary Response File grid. 
The algorithm is described in more detail in lButler. et al.l (l2005al ). We model the LRF versus 
PI energy bin as a Gaussian with an energy dependent width, a^E) = 29.4(i?/[lkeV])°'^^^ 
eV. This functional form provides an excellent fit to the core of the XRT LRF. We estimate 
the significance of each line set by applying the Likelihood Ratio Test (LRT) (also referred 
to as the Ax^ test) to the difference in observed x^ values prior to and after adding the 
emission line component. These significance estimates are then checked by comparing to 
the distribution of Ax^ produced by a Monte Carlo (MC) simulation (Section [3]). We quote 
significance values in units of normal distribution cr's (i.e., a null hypothesis probability of 
2.7 X 10-3 is 3a, 5.7 x 10"^ is 5a, etc.). 

For a spectrum containing ~ 500 counts, the typical energy binning oversamples the 
LRF full-width at half max by a factor ~ 2. By fixing the number of counts in each spectra, 
it is a simple task to determine MC line significances for the sample as a whole. Such a 
determination for each spectrum would be prohibitively time consuming. To be sure that we 
are not missing features due to the 500 count constraint, we have also dyadically grouped 
(2,4,8, spectra together, etc., up to the full integration) and searched the spectra for each 
burst. We impose a tighter trigger criteria (see Section 14. 4p on these data. In section 15. H 
we discuss a line search for the integrated spectra of bright (> 10 cts/s) XRT flares. 

Figure [3] displays our triggers for 1-5 emission lines from the blind search through the 
full data set. Searching 1152 spectra, the trials probability for finding one spectrum with a 
4.5a trigger is 99% confidence. Given the sample median, lower and upper quartile values for 
the observed single-trial trigger significances found — 1.39a, 0.94a, and 1.83a, respectively, 
4.5a is also 3 interq uartile lengths from the upper quartile, which is a common measure for 
strong outliers (e.g., lDevorelll995l ). We find a handful of triggers near or beyond this value, 
which occur in the spectra of 3 bursts (e.g.. Figure [3]). As mentioned above, the significance 
estimates used here come from the LRT test. We also check the validity of this procedure by 
simulating 10^ spectra with 500 counts each and searching for lines in the simulated spectra. 
For the simulation, we take F = 2 and Nh = 10^^ cm"^ as representative values for the 
entire XRT sample. The resulting (false) trigger rate for the most significant subset of 1-5 
emission lines from the search is plotted in blue in Figure [31 There is a clear departure at 
low significance from the model (dashed curve) for the LRT. Importantly, the departure is 
negligible and conservative for si gnificances >, 3'7. The validity and l imitations of the LRT 
are discussed in detail elsewhere ( lYaqooblll998uProtassov et al.ll2002l ). 
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Fig. 3. — The number of triggers (red lines) for 1-5 emission lines in XRT spectra containing ~ 500 
counts versus their estimated significance. A branch in the XRT curve is seen depending on whether 
or not 050714B, 050822, and 060218 are included in the sample. The bottom red curve does not 
include these events and is roughly consistent (KS-test significance ~ 3a) with the Monte Carlo 
false event rate (blue line). Ignoring the data from the anomalously bright GRB 060218 — which 
comprise ~ 30% of the spectra under study — the excluded 2 bursts correspond to 4% of the the 
full sample spectra. The Monte Carlo curve is determined from 10^ simulations of a fiducial 500 
count spectrum (Section [3]). The dotted black line is the expected false event for the LRT (i.e., the 
cumulative normal distribution). 
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Fig. 4. — The fraction (= efficiency) of 10^ simulated spectra with lines of varying signal to noise 
(S/N) where the trigger algorithm found the lines at their true locations and with consistent S/N 
values. The efficiency is near unity for 1-5 emission lines, each with S/N ^ 3. The efficiency drops 
sharply for multiple lines with S/N ^ 2. 



We also perform simulations to determine the efficiency of the algorithm for detecting 
lines known to be present. Starting with the same continuum model used for the false trigger 
tests, we add lines with S/N = 3 (comparable to low S/N values from lines detected previ- 
ously) randomly placed in energy over the 0.3-5.0 keV band. Figure H] displays the efficiency 
(fraction of iterations where the lines are detected at their true locations) of the algorithm 
as the assumed line S/N is varied in the MC The efficiency drops sharply for S/N ^ 2 but 
remains near unity f or S/N ^ 3. We note that the trigger software autonomou sly finds (see, 
Butler, et al.l l2005af) each of the 5 em ission lines claimed for GRBs 011211 ( Reeves et al. 



2OO2I ) and 030227 (jWatson et al.ll2003l ). There is no clear dependence of the efficiency on the 
line energy, however, the restricted energy band places clear limits on the chemical species 
to which the search is sensitive. For nearby events {z <^ 0.4), the upper energy cutoff of 5.0 
keV results in non-detection of the Fe- group elements. The low-energy cutoff can inhibit 
detection of light element lines (e.g., Ne at 2; ;^ 2, Si at 2 ;^ 5). From the simulations, we 
estimate that LRT significances in a units are accurate to ~ 20% for S/N ^ 3 lines with 
o'hRT ^ 3. Potential unresolved emission lines from astrophysically abundant species are 
detected with high efficiency over a broad range of redshifts. In principle, the search is also 
sensitive to broad emission lines, which can be built up from the superposition of multiple 
narrow lines. We now discuss the individual triggers in reverse chronological order. 
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4. Line Emission Detections in Soft Portions of Three X-ray Afterglows 

4.1. GRB 060218 
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Fig. 5. — (A) The X-ray light curve for the GRB 060218 afterglow. "Hardness" is the difference 
over the sum of the 1.5-8.0 keV counts and the 0.5-1.5 keV counts. (B) Power-law fits to the 
afterglow. 



The flux from this event was extraordinarily high. The light curve displays a prominent 
hard to soft evolution during a broad temporal rise, followed by a rapid decline and leveling 
off (Figure E]). The soft portion of the afterglow (after t ~ 1 ksec) accounts for much of the 
tai l of the F ~ 2 — 3 distri bution in Figure [21 A preliminary spectral analysis is reported 
in ICusumano et al.l (l2006bl ). We group the WT (PC) mode counts into 355 (2) spectra, 
each containing ~ 500 cts. Due to the large data mass, we present our fit results in plot 
form (Figure 03) for these finely resolved data, providing tables below for the spectra with 
~ 16,000 cts. 

Overall, the finely time-sliced data are well fit {x^ /v = 12149.5/12291 for 356 spectra) 
with absorbed power-law's, with Nh = (3.8 ±0.4) x 10^^ cm~^ on average (consisten tly larger 
than the Galactic value A// Galactic = l-H x 10^^ cm~^; iDickey fc Lockmanlll990l ). Toward 
the end of the WT mode data at (t = 2.317 — 2.326 ksec), the power-law fit becomes poor 
(x^/u = 52.59/32, with F = 2.3 ± 0.2, Nh = (3.9 ± 0.6) x lO^^ cm'^, and unabsorbed flux 



/ = (3.1 ± 0.3) X 10~^ erg cm~^ s^^ [0.5-10 keV]). The fit is improved at 4a significance 
(Ax^ = 24.69, for 4 additional degrees of freedom) with the addition of 2 emission lines 
(Figure [7]). The best fit lirie ene rgies are 0.82±0.04 and 1.01 ± 0.04 keV. For redshifts near 
z = 0.033 (JMirabal et al.ll2006l ). the lines can be identified with K-shell transitions in H 
and He-like species of Ne or L-shell transitions from Fe-group elements. The equivalent 
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widths and luminosities are 205±42 and 145± 30 eV and 3.6 ± 0.2 and (1.9 ± 0.2) x 10"^° 
erg cm~^ s~^, respectively. The blackbody model provides a poorer fit to the data in this 
{yd /^ = 82.96/32) and the other time regions. A power-law plus blackbody model provides 
a mediocre fit {x^/u = 42.94/30). 

As we discuss below, the integrated spec trum is also quite well fit by a thermal MEKAL 
(JMewe. Gronenschild. fc van den Oordlll985l ) plasma model {x^/v = 26.90/30) with kT = 
0.44 ± 0.06 and normalization (see note in Table 3) 3.0 ± 1.3 for solar abundances. The 1-a 
lower limit on the abundance is 0.2 solar, with an undefined upper limit. The line emission 
is dominated by L-shell Fe, with an Fe to Ne abundance ratio 1.5^*^0^7^ 5 relative to solar. 



— 10 



47 



w 

'w 
o 



10 
10 
10 
10 
10 



46 
45 



> 0.2 
H 0.1 



44 



43 



42 



3 h 



■+H- 



-I— 1--1- 



PL Model — I- 

BB Model — ^ 

Emission Lines — *- 



-i I I I I I L 



J_ 



10"" 
Time Since Burst [s] 



-i I I I I I L. 



10^ 



t3.0±0.3 +-0.6±0.2 



Radius 



....J H 10 



10-" 



- 10 



12 



^11 



10^ 



Fig. 6. — A spectral deconvolution for GRB 060218 into soft and hard emission components. 
The plotted power-law model components (Flux and photon index F) are for the power-law plus 
blackbody model. The blackbody radius is measured in cm. 
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WT Mode Spectrum with 2 Emission Lines 




Energy [keV] 

Fig. 7.— The XRT WT Mode spectrum from 2.317 to 2.326 ksec after GRB 060218. Two 
emission lines at 0.82 and 1.01 ke V (green curve), possibly associated with L-shell transitions in Fe 
at z = 0.033 (IMirabal et al.ll2006l ). improve the power-law fit (red curve) at 4(T significance. In the 
residuals panel, the green curve shows the way the power-law residuals (black curve) are fit by the 
power-law plus lines model. 



There is strong evidence that the soft emission component is present throughout an 
extended portion of the observation. Only ~ 50% of the dyadically grouped spectra con- 
taining > 500 cts are found to be adequately fit by the simple absorbed power-law models. 
Particularly for the spectra containing 16,000 or more counts, a soft component near 1 keV 
is strongly required (power-law model alone rejected at > 90% confidence in 11 of 12 spectra; 
Table 1). With the inclusion of the blackbody component or emission lines in addition to 
the power-law, each fit to the ~16,000 cts spectra yields x^/^ ^ 1- 

Figure M shows the results of the spectral deconvolution of these data into power-law 
and soft components. The power-law component displays a smooth hard to soft evolution in 
time. In the case of the emission line model, the lines in the 0.8-1.0 keV energy range, are 
detected in multiple time intervals (Table 2), with high (> 5cr) significance triggers coming 
after t ~ 1 ksec. In addition to the lines at E <^ 1 keV discussed above, triggers in the 3-4 
keV range are found (Table 2), possibly associated with H- or He-like Ar and Ca. 

The soft-component luminosity L in both cases appears to peak after the power-law com- 
ponent. The blackbody temperature declines mildly after t = 580s as (0.19±0.01)(t/580s)*-"°'^^^°'°^'' 
keV to kT = 0.105 ±0.004 keV at t = 1.52 ksec, then remains constant and possibly declines 
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at late time. Using the standard formula for the blackbody radius: 



R 



BB 



2.78 



lO^^erg s 



-1 



1/2 



kT 



0.1 keV 



X lO'^cm, 



and z = 0.033 (JMirabal et al.ll2006l ). the blackbody radius increases, peaks, then decreases 
with the time dependences shown in Figured Here and throughout, we assume a cosmology 
with (h,r2m,r2A) = (0.65,0.3,0.7). We find consistent behavior independent of whether or not 
we fix Nfj or kT in the fi ts. Consistent values for the initial blackbody radius are reported 
by (jCampana et al.ll2006l ). However, they find that the radius continues to grow in time 
during the XRT observation without peaking, a scenario that we find to be ruled out at the 
~ 7a level. 
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Fig. 8. — Power-law models (red curves) fail to describe the soft emission in the integrated WT 
mode (t = 160s-2.8ksec) and PC mode {t = 5.9ksec-3.1Msec) spectra for GRB 060218 (Table 3). 
(A) Top two plots. Both data sets are better fit with the addition of 4-5 emission lines (green 
curves, also plotted for the power-law model residuals): WT mode, 4 lines, 20. Go", Ax^ = 460.03, 
for 8 additional degrees of freedom; PC mode, 5 lines, G.9a, Ax^ = 75.37, for 10 additional degrees 
of freedom. (B) Bottom two plots. The fits to both data sets are also markedly improved with the 
addition of a blackbody continuum component (Table 3). 



Figure [8] shows the integrated WT and PC mode spectrum for the afterglow. Continuum 
spectral fits are reported in Table 3. Prominent residuals are present in both spectra near 1 
keV. The WT and PC mode data can be fit with 4 emission lines at energies 0.80 ± 0.03, 
0.9 ± 0.03, 1.00 ± 0.03, and 1.1 ± 0.03 keV in addition to the power-law continuum. The 
WT mode trigger is the most significant trigger in the entire sample by a large margin (see 
Section W^ . The line luminosities are ^ 2 x 10^^^ erg cm~^ s~^ in the WT mode spectrum 
and ~ 10^ times fainter in the PC mode data. Oppositely, the equivalent widths increase by 
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a factor ~ 5 from ~ 20 eV in the WT mode data, indicating that the hne fluxes decrease less 
slowly than the continuum flux (see, also. Figure [T6|l . If we relax the upper limit on the line 
energy search, there is a weak (~ 2a) trigger on a possible K-shell Fe- series line at -E = 6±1 
keV in the PC mode spectrum, with an equivalent width of order 1 keV (e.g., residuals in 
Figure [8j3). The further addition of lines to either spectrum does not significantly improve 
the fits. 

The poor apparent x'^/'^ = 1280.53/760 in the case of the WT mode data probably 
should not be taken too seriously, because the flux is high and the inclusion of a small ^ 5% 
systematic error component leaves the fit statistically acceptable. The spectra can also be fit 
with a thermal plasma component (assuming solar abundances) in addition to the power-law 
(Table 3)-which reinforces the notion that some of the soft emission may be in the form of 
discrete lines-or with a blackbody in addition to the power-law. 

Comparing the hypothesis of emission lines to that of a blackbody (in additional to a 
non-thermal continuum) for GRB 060218, we find that the event is a particularly strong line 
candidate. The power-law plus two line model is an excellent fit to the 500 cts spectrum dis- 
cussed above (x^/^ = 27.9/28), whereas the power-law plus blackbody model fit is mediocre 
(x'^/u = 42.94/30.) 



4.2. GRB 050822 

This ev ent is a clea r outh er in Figure [2l A preliminary spectral analysis has been 



reported by iGodet et al.l (120051 ) . The spectrum displays a gradual softening in time which 
increases abruptly with the impulsive flare at t ~ 430s (Figure (HP). The continuum spectrum 
during the flare (430s to 550s after the burst) is well- fit by a blackbody (Table 4), with a 
decreasing temperature in time kT = 0.26—0.17 keV (observer frame). There is an indication 
that the blackbody radius from Equation [1] may increase in time by ~ 50%, starting at 
i?BB = (2.7 ± 0.5) X 10^2 cm at t = 430s and peaking at Rbb = (4.5 ± 1.3) x 10^^ ^.^ ^t 
t = 490s, assuming z = 1.2 (see below). 

A power-law model also provides an acceptable fit to the data during the flare, except 
during the penultimate time interval, 489.5 to 509.4s (interval ^14 in Table 4; x^/z/=54. 71/33, 
rejectable at 99% confidence). The blackbody model fit is marginally better in interval 
#14 (x^/i/=51. 06/33, rejectable at 98% confidence) than the power-law fit. An added 
set of 5 emission lines improves the power-law fit during interval #14 at 4.4cr significance 
(Ax^ = 40.56, z/ = 10; Table 4). The same 5 lines improve the blackbody fit at 2.8cr signif- 
icance (Ax^ = 25.03, z/ = 10). The continuum fit residuals are narrow (Figure fTUj) : neither 
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the power- law nor the blackbody fit for interval #14 is improved by adding an additional 
blackbody or power-law continuum component. 
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Fig. 9. — (A) The X-ray hght curve for the GRB 050822 afterglow. There is a strong flare at 
t ~ 500s, during which the spectrum becomes very soft. (B) Spectral evolution for GRB 050822. 
A dramatic increase in the spectral curvature near 400s is clear from a jump in the photon index 
and best-fit absorption column. 



We test for the presence of the line set found for interval #14 by fitting the best the line 
set model to the data from the other time intervals. We allow one fioating normalization, 
with the other line parameters fixed relative to this normalization. If the line-set from interval 
#14 is present during the fiare outside of interval #14 (i.e., in intervals 11, 12, 13, and 15), 
the combined fiux from the 5-line set must be lower by a factor > 4 (90% confidence). Thus, 
the overall or relative fiux of each line appears to be changing in time. If we take the ~ 20s 
period during which the line triggers are found to be most sigificant as representative of the 
period of active line emission, then the possible detection appears very unusual and probably 
not credible given the models for line emission (Section 15.31 below). We note that there is a 
weak trigger on the line set at earlier times (intervals #3-#5), with a mean fiux consistent 
with that found for interval #14 but also consistent with zero at the l-a level. 

By hand, we combine the PC and WT mode spectra in the tail of the fiare (see Figure 
[TOl) and use the PC mode redistribution matrix to perform a fit. The data in this interval 
(#19 in Tables 4 and 5) are also poorly fit by either continuum model and show modest 
significance evidence for 5 emission lines. Thus, the very short duration of the line trigger 
may be due to the detection threshold and may not refiect the physical emission timescale. 

The redshift is currently unknown for this event. The line centroids during interval 
#14 (Table 4) are too closely spaced to allow for an association with H-like ions from light 
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metals only (e.g., iReeves et al.ll2002l : IWatson et al.ll2003l ). However, if we allow the He-like 
species, identifications are possible and non-unique. One possible association for the 5 lines 
in Figure crop is: Ar XVII, S XVI, S XV, Si XIV, and Si XIII at z = 1.2. From Table 5, the 
lines at E = 0.81, 0.91, 1.04, and 1.23 keV are detected in multiple time intervals. A line 
near 3.5 keV is detected in two time intervals and could be associated with H-like Co or Ni 
at z = 1.2. A power-law plus MEKAL model provides a poor fit due to the closely spaced 
emission lines. 
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Fig. 10.— (A) In interval #14 during the GRB 050822 flare, emission lines at 0.81, 0.91, 1.04, 1.23, 
and 1.49 keV (green curve), possibly associated with Si XIII, Si XIV, S XVI, S XV, and Ar XVII 
at z = 1.2, improve the power-law model fit (red curve) at 4.4(7 significance (Table 5). (B) The 
same but for a blackbody continuum model. The lines are 2.8o" significant. For the residuals panel 
in each plot, the green curves show the way the continuum model residuals (black curves) are fit 
by the continuum plus lines model. 

There is one minor caveat relevant to our analysis of this event. We include two detec- 
tor columns (RAW X columns #290 and #291) in the WT mode data which are normally 
discarded in the standard xrtpipeline processing. These are neighboring columns (sec- 
tions of which have been found to behave anomalously) to columns which were damaged as 
the result of a possible micro-meteorite impaccl. We find that the continuum spectral fit 
parameters are consistent whether the columns are retained or rejected. However, because 
the two columns are near the center of the source extraction region for a large part of the 
observation, the loss in source flux associated with rejecting the columns mildly reduces the 
line trigger signiflcance to 3.5a. Here the time extraction is broadened to 480-512s, in order 



'http://swift.gsfc.nasa.gov/docs/swift/analysis/xrt_digest.html 
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to get 510 counts in the spectrum. 

Comparing the hypothesis of emission hues to that of a blackbody (in additional to a 
non-thermal continuum) for GRB 050822, we find that this event, like GRB 060218, is a 
strong line candidate. The power-law plus five line model provides an excellent fit, {x^ /v = 
14.15/23), whereas a power-law plus blackbody improves the power-law fit very little. We 
note that the line sigificance degrades (from 4.4o" to 2.80") if the underlying continuum is 
assumed to be a blackbody. 



4.3. GRB 050714B 
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Fig. 11. — (A) The X-ray light curve and hardness plot for GRB 050714B. (B) There is strong 
spectral evolution during the early X-ray afterglow. 



Also a clear outlier in Figure [2l the WT mode spectrum for the time interval 157.4 
to 189.7s (#1 in Table 6) corresponds to a rapid decline in the X-ray light curve and is 
well-fit by a blackbody. Alternativel y, the spectrum c an be modeled by a power-law with 
photon index F = 5.9 ±0.4 (see, also. lPage et al.ll2005al ). The spectrum during the flare (PC 
mode; t ^ 275 — 525s) is poorly fit by either a blackbody {x^/i' = 97.19/35) or a power-law 
ix^/^ = 76.78/35). The PC mode spectrum from the start of the flare until 58.3 ksec, can 
be marginally well fit (x^/i^ = 48.64/33) by a power-law plus blackbody. In Table 6, we 
also extract the PC mode data during the flare only (interval #3) and demonstrate that the 
same result holds there, although with larger error bars. Between regions #1 and #3, the 
blackbody temperature in the observer frame decreases from 0.2 to 0.1 keV, and the radius 
increases from i?BB = (2.7 ± 0.5) x 10^^ cm to (1.2 ± 0.7) x 10^^ cm (Equation [1]), assuming 
z = 2.66 (see below). 
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Alternatively, the power-law fit to interval #2 is improved at 4.2cr significance with the 
inclusion of 4 emission lines {x^/v = 41.69/27; Ax^ = 35.09, z/ = 8; Figure [T2l) . The fit 
to region #3 is improved less with the addition of lines, however there is a small number 
of counts in that hand-selected region, and the line locations are consistent with those from 
region #2 (Table 7). There are a number of plausible line identifications for region 7^2. For 
the lines at 0.91 ± 0.03, 0.76 ± 0.05, 0.56 ± 0.05, and 1.12 ± 0.05 keV, one association is 
made with H-like species Ar XVIII, S XVI, Si XIV, Ca XX, respectively, at z = 2.66. A 
thermal (MEKAL) plasma with normalization (see. Table 3) 0.1 ±0.05 for solar abundances 
and kT = 1.0 ±0.2 keV at a similar redshift [z = 2.4), in addition to the power-law, provides 
an acceptable fit to the data ix^ /v = 44.48/33). 




Energy [keV] 
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Fig. 12.— (A) During tfie GRB 050714B flare at t ~ 400s, emission lines at 0.56, 0.76, 0.91, and 
1.12 keV (green curve), possibly associated with Si XIV, S XVI, Ar XVIII, and Ca XX at z = 2.66, 
improve the power-law model fit (red curve) at 4.2a significance (Table 7). (B) The power-law fit 
is also improved significantly with the inclusion of a blackbody (Table 7). The residuals panels in 
each plot show how the additional model components (green curves) fit the power-law continuum 
model residuals (black curves). 



Comparing the hypothesis of emission lines to that of a blackbody (in addition to a 
non-thermal continuum) for GRB 050714B, we find that this event is a compelling line 
emission candidate. The power-law plus line model yields an excellent fit {x^ /v = 41.69/27). 
However, this is little improvement over the simpler power-law plus blackbody model {x'^/iy = 
48.64/33). The better fit of the MEKAL plasma model {x^/u = 44.48/33) for the same 
number of degrees of freedom appears to reinforce the line emission hypothesis. 
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4.4. Detections in the Dyadically Grouped, > 500 Count Spectra 

The sample of spectra under study is roughly doubled when we consider the additional 
spectra formed by dyadically grouping the ~ 500 count spectra into more fully-exposed 
spectra containing > 500 counts. We observe that the inferred line significance is clearly a 
function of the number of counts within the spectrum, with the more fully exposed spectra 
showing a higher probability of false trigger. The median significance (from the LRT) for 
1-5 emission lines for the full sample is 1.71a. The lower and upper quartiles are 1.21cr and 
2.23(j, respectively, implying a limit of 5.3o" for a strong outlier. We find that only 3 bursts 
exhibit spectra with line triggers at the > 5a level: GRBs 060218, 060202, 050822. The 
maximal significances are 20.6cr (t = 160s-2.8 ksec; 27,620 cts), 6.5o" {t = 150s-l ksec; 27,620 
cts), and 5.2a [t = 482.9 - 525.5s; 1050 cts), respectively GRBs 060218 and 050822 are 
discussed in detail above. Here, we briefly present the GRB 060202 spectrum. 

The WT mode data during the first 1 ksec of GRB 060202 show two temporal de- 
clines separated by a plateau, with mild spectral evolution (Figure [T^ then and after- 
ward. The spectrum in the 0.6-10 keV band from 150s to Iksec after the burst is un- 
acceptably fit by a blackbody model ix^ ju = 2511.64/507). Less poor is the absorbed 
power-law model fit {x^/v = 644.54/507), with the following best-fit parameters: Nh = 
(5.0 ± 0.1) X 10^^ cm~^ (significa ntly in excess of the Galactic column: A^_h Galactic = 5 x 10^° 



cm 



Dickev fc Lockmanlll990h . / = (2.04 ± 0.02 ) x 10"^ erg cm"^ s"^ [0.5-10 keV], and 



r = 2.32 ± 0.02. Consistent values are reported by iMorris et al.l (120061 ). The fit is improved 
at 6.5(7 significance with the inclusion of 5 emission lines (Figure [TH Ax^ = 69.06, for 10 
additional degrees of freedom). The best fit line energies are: 0.9 ±0.03, 1.01 ±0.03, 1.12 ± 
0.07, 4.7 ±0.07, and (1.19 ±0.04) keV. The line equivalent widths are 17 ± 3, 17 ± 2, 17 ± 
2, 50 ± 15, and 10 ± 2 eV, respectively. The line luminosities are 1.23 ± 0.01, 1.4 ± 0.1, 
1.6 ± 0.1, 0.5 ± 0.2, and (0.9 ±0.1) x 10"^^ erg cm'^ s'^ respectively 

From Keck optical spectroscopy, we have determined that the host galaxy redshift is 
z = 0.783. The X-ray lines can be indentified with Al XIII, Si XIII, Si XIV, Ni XXVIII, and 
S XV, respectively. In an initial draft of this paper, we had guessed z = 0.4 for this event. 
Altough the prediction was not borne out in detail, the complex of X-ray lines near 1 keV 
did correctly indicate a moderately low z. No lines are detected after 1 ksec at the > 2.5cr 
significance level. At the 3a level, the line luminosities are below 1.6 x 10^^^ erg cm~^ s~^ 
and the equivalent widths are below 300 eV. 

Alternatively, we can model the soft component by adding a blackbody with kT = 
0.18 ± 0.02 keV (observer frame) and luminosity (7.4 ± 2.4) x 10^^° erg cm~^ s~^, implying 
a radius (3.0 ± 0.8) x 10^^ cm (Equation [T]). The reduced x^ is marginally lower than for the 
power-law plus lines model {x^ /v = 570.68/505). The absorption is Nh = (5.9 ± 0.3) x 10^^ 
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cm~^, and the power-law model parameters are / = (1.8 ± 0.1) x 10~^ erg cm~^ s^^ [0.5-10 
keV], and F = 2.20 ± 0.04. The thermal component can also be modeled by a MEKAL 
plasma. However, the implied abundances are zero, indicating that the data prefer an 
additional continuum (rather than discrete) component. 
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Fig. 13.— (A) The X-ray light curve and hardness plot for GRB 060202. (B) The absorbed 
power-law model fits. 
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Fig. 14. — (A) The absorbed power-law model fit (red curve) to the WT mode data 150 to 1000s 
after GRB 060202 is improved significantly with the addition of 5 emission lines (green curve) at 
0.9, 1.01, 1.12, 1.19, and 4.7 keV, possibly associated with Al XIII, Si XIII, Si XIV, Ni XXVIII, 
and S XV at z = 0.783. (B) A slightly better fit is found by adding a blackbody to the powerlaw. 
The residuals panels in each plot show how the additional model components (green curves) fit the 
power-law continuum model residuals (black curves). 
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Comparing the hypothesis of emission hnes to that of a blackbody (in additional to 
a non-thermal continuum) for GRB 060602, we find that the blackbody model provides a 
better fit {x^ /v = 570.68/505) than the more complex power-law plus lines model {x^ /v = 
575.45/497). 



Discussion 



5.1. The Spectra of Other Bright Flares 

Given the detection of soft X-ray components in two bright flares, it is possible that a 
correlation may exist between this emission and the X-ray flaring. To explore further this 
possibility, we extract and fit spectra by hand to the integrated counts of 19 bright (peak 
count rate > 10 cts/s) flares (Table 8). Because a blackbody model provides a poor flt to 
each spectrum, we present only the power-law model flts (Table 8). The F values for GRBs 
050714B and 050228 stand out from the other values, most of which are F ~ 2. (There is 
also one hard outlier, GRB 050820, with F = 0.99 ± 0.05.) We therefore conclude that these 
very soft spectra are uncommon in the flares. 

A few of the flares, in addition to those discussed above, exhibit marginal significance 
(> 3a) evidence for lines on top of their non-thermal spectra (Table 9): GRBs 050502b, 
051117A (during two flares) , 060124 (flrst flare, with weak detections in the next). For 
GRB 060124 at z = 2.296 JProchaska et al.ll2006h . there is a possible Fe XXVI line near 
1.9 keV and several lines at higher energies possibly associated with recombination of Fe- 
group elements. Although no redshift is known for GRB 050502B, the lines may also be 
associated with the Fe group elements: 1.44 keV— Fe XXVI, 1.69 keV— Ni XXVII, 4.48 
keV — ?, 2.14 keV — Fe ionization? One possible association for the lines in GRB 051117A at 
energies ~ 1.8, 2.8, 0.6, 1.4, 0.9, and 1.1 keV, is with H- or He-like Ca, Fe, Ar, Ne, Si, and 
S, respectively, at z = 1.2. We expect ~ 1 trigger at ;^ 3cr in 22 trials. The actual number 
of triggers (6) suggests a modest increase in the line emission probability during the bright 
flaring episodes relative to the quiescent periods (e.g.. Figure [3]). 



5.2. Blackbody Emission: Shock Breakout Through the Progenitor Star? 

The optical/IR emission from GR B 060218 has yielded unequivocal evidence for an un 



deriving type-Ic supernova-SN 2006aj (IMirabal et al.ll2006l: ISqllerman et al.ll2006l : iModjaz et al. 



20061 ). It has been suggested recently (ICampana et al.ll2006l ) that the unusual early X-ray 
emission is due to the propagation of a radiation dominated shock through a wind or H en- 
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velope surrounding the progenitor star. In this picture, the early X-ray hght curve is broad, 
then rapidly declining, due to light travel time effects across an aspherical shock shell of 
thickness ~ R/Tc ~ 300s. 



The shock can be due to the GRB itself (ITan. Matzner. fc McKeell200ll : IColgatelll974[ ) 
or due to mildly relativistic material expected to form and surround t he GRB jet as it 
punches through the He core of the prog enitor star at Rue ^ 10^^ cm ( JMeszaros fc Rees 
200ll : iRamirez-Ruiz. Celotti. fc Reed 120021 ). After the shock punches through an optically 
thick layer around the He core and interacts with the envelope, it leaves behind an expanding 
photosphere. The radius at which the light is abl e to escape the photosphere aro und a carbon- 
rich, Wolf-Rayet star is i?p = 1.5 x lO^^^ cm jTan. Matzner. fc McKe3l200l[ ). The scatter 

in observed mass loss ra tes (M_4 in units of 10~^Mp) yr~^) and wind velocities [v^^g in units 
of 10^ cm/s) is ~ 1 dex (JKoesterke fc Hamannlll995l ) and could lead to a ~ ±1 dex range in 
Rp. Due to the energy input (~ 10^^ erg) the expected photospheric blackbody temperature 
is ~ 0.3 keV. 

In the course of our line search, we find spectra from 3 other events which are soft 
and well fit by mode l s containing a blackbody. Aside from no clear detection in the UV 
( IGronwall et al.ll2005l : iPage et al.ll2005bl : iBlustin et al.ll2006l ). these detections are similar to 
GRB 060218 in a number of ways. In each case the observation epoch is ~ 1 ksec, and the 
blackbody radii are Rbb ~ 10^^ — 10^'^ cm. We also find similar temperatures (~ 0.1 — 0.5 
keV) in each case, which appear to decrease in time mildly. In all cases, the emission is 
transient. The Rbb/c are similar to the light travel time from the GRB to the observation 
epoch, which would imply a simultaneous GRB and SN in either scenario discussed above. 
The radii are at least 2 orders of magnitude smaller than the distance expected for the highly 
relativistic (F ~ 10) external shock at the observation epoch, ruling out an association with 
the external shock. 
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Fig. 15. — The isotropic X-ray luminosity at t = 10 hours in the host frame versus the isotropic 
prompt energy release in 7— ray's. Fluence data for Swift (typically in the 15-150 keV band) 
are taken from the table at http://swift.gsfc.nasa.gov/docs/swift/archive/grb_table, and these are 
translated to -E^^iso by multiplying by 47rL)£/(l + z) only. The bursts studied here are shown 
encircled in blue. The parallel track set by the Swift bursts relative to bursts prior to Swift likely 
reflects only the narrow BAT band pass and our non-application of a k-correction, which apparently 
leads to systematic underestimates of -E^^iso- Note the shift in ii^-^^iso for GRB 06 0218 depending 



on wh ether we use the tabulated value or the more carefully calculated value from lCampana et al 
(|200d ). 



An additional clue to the mechanism which produces the soft component may come 
from a simple comparison of the GRB and shock kinetic energies. We can infer the shock 
kinetic energy using the X-ray afterglow luminosity at t = 10 hours in the host frame, which 
is expected to scale linearly with the shock energy due to an indepe ndence on the density and 
a mil d dependence on the shock microphysical parameters (e.g., iBerger. Kulkarni. fc Frail 
2003h . GRB 060218, like GRB 980425/SN 1998bw and GRB 031203/SN 20031w, is found to 
track other GRBs from Swift and previous missions (Figure [TSll . implying a roughly similar 
efficiency for the conversion of shock kinetic energy into 7— ray's, but with a substantial 
sub-luminosity evident in both quantities. The other bursts discussed above also follow the 
trend Lx,iso oc -E^,iso- 



On the other hand, using the fact that the duration of the photospheric expansion phase 
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appears to be ~ 10^ — lO^s for each burst, we can estimate the ratio of blackbody fluence 
to prompt GRB fluence (Table 10). The prompt 7— ray energ y release for GRBs 060202, 
050822, and 050714B is S^ = 2.4 ± 0. 2 jHuUinger et al.lbooeJl . 3.4 ± 0.3 jHulhnger et al 



2005|), and (0.7±0.r 



X 10 ^ erg cm ^ (JTueller et al 



20051 ). respectively, in the 15-150 keV 



band. The ratios are found to be: 0.3 ± 0.1, 0.03 ± 0.01, and 0.2 ± 0.1, respectively. The 
ratio value for GRB 060218 is 4 ± 1. The soft fluence clearly correlates with the prompt 
7— ray and shock emission. The large dispersion over more than two orders of magnitude 
could be explained by the diversity of observed progenitor wind outflows and the strong 



dependence of the available en ergy E on the mass loss history: E 
fiTan. Matzner. fc McKe"3l200lh . 



r^R 



Mc^ 



OC 



'M_\2 



Because wide variations in the thermal energy release apparently do not alter the Lx^iso oc 
Ey^iso relation, we conclude that the soft component is more likely to have been produced 
by the cocoon surrounding the GRB jet as it escape the progenitor He core than by the 
GRB jet itself. A radiation dominated cocoon is not expected to appreciably affe ct the 
standard internal shocks or afterglow emi ssion (IRamirez-Ruiz. Celotti. fc Reed l2002l ) . On 
the issue of the faintness of GRB 060218, iFan. Piran. fc Xd ([20061) discuss the possibility 
that GRB 060218 was an almost "failed GRB," with a large fraction of the GRB energy not 
making it out of the progenitor star and winding up as thermal photons in the cocoon. 

One additional insight which comes from the extreme proximity of GRB 060218, it's 
sub-luminous 7— ray emission, and the relative faintness of the soft X-ray components in 
3 other bursts is that this emission could be relatively common. Bursts like GRB 060218, 
where the spectrum is dominated by the soft component could outnumber normal GRBs 100 
to 1. It may be a selection effect that the blackbody radiation happens to peak where the 
detector response is also maximal at -E ~ 5kT ~ 0.5 — 2 keV for each of these events; soft 
or harder spectra due to breakout at larger or smaller radii, respectively, may also occur. 



5.3. Constraints on the Line Emission Mechanisms 



In the previous section we discussed the soft X-ray emission in terms of the blackbody 
flts. For 3 of the 4 bursts in Section IH the data appear to be flt better with an emission 
line model. Clear line associations are possible for GRB 060218 and possible association 
are presented for the other bursts lacking redshift measurement. For two of the bursts, 
the line component is well flt by a thermal plasma model containing only two parameters 
(temperature and normalization). Because, the continuum under the line s is domi i iated by a 
power-law, the emitting plasmas may well not be in thermal equilibrium. iLazzatil (120031 ) has 
discussed the relative unlikelihood that GRB X-ray lines could be due to thermal plasmas. 
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For simplicity in fitting, however, the MEKAL code generates approximately correct line 
locations and sensible flux ratios for a broad range of astrophysically abundant elements. 
The code has beeii shown to produce similar lines as from photo-ionization models (e.g., 
Watson et al.ll2003l ). at the level of spectral resolution appropriate to Swift XRT data. In 
the case of GRB 050822, the MEKAL model cannot produce the closely spaced lines, and this 
may argue for a non-equilibrium ionization state, for two plasmas with differing ionization 
states, or for plasmas moving at different speeds (^ 0.1c) relative to the observer. 

In none of the cases do we detect lone, statistically significant emission lines. Rather, 
we find sets of lines, typically fairly closely spaced in energy, which we associate with ionized 
light metals or with L-shell transitions in Fe (GRB 060218). The line sets we find are 



perhaps most similar to the lig ht metal lines ( Reeves et al.ll2002l : IWatson et al.l l2003l ) and 
soft, possibly thermal excesses (IWatson et al.l 120021) found in XM M data. A detection of 
such lines has been claimed for Chandra as well ( JButler et al.ll2003l ). 



Figure [T6l adapted from lButler. et al.l (J2005bl ). displays the hne triggers in the context of 
lines claimed from previous missions. For comparison and for estimation of upper limits for 
the ~ 90% of spectra not showing significant evidence for hues. Figure [T71 displays equivalent 
width values versus time for marginal detections in the full dataset of Swift bursts with 
measured redshift. It is difficult to use the triggers or upper limits to rule on the reality of 
the historic lines, because there are few well-exposed XRT spectra at t ;^ 0.1 day. However, 
observatio ns with Chandra pl ace stringent limits on the late-time lines (down arrows in 
Figure flGj iButler. et al.ll2005bl ). The upper limits from the non-detections {EW ^ 1 keV at 
t < 0.1 day) in the Swift XRT sample are roughly consistent with the detections. 

The most striking feature of the triggers is that the possible lines are emitted at dra- 
matically earlier times (t < 0.1 day) than in previous cases. Prior to Swift, it was suspected 
that GRB lines w ould be very difficult to forrn at ea rly times (aside from exotic possibilities 
suggested by, e.g., iGou. Meszaros. fc Kallmanll2004l ) due to the overwhelming flux from the 
bright X-ray afterglow. The solid lines in Figure [16] show the rapid decrease expected at 
early t imes for two photo-ionization models (JBallantyne fc Ramirez-Ruizll200ll : iButler. et al. 



2005bl ). The data from GRB 060218 may roughly follow this expectation, but the other data 
likely do not. 

One of the big surprises from Swift was the departure of the early X-ray afterglow light 
curves from late time extrapolations (early flaring, anomalo usly flat and rapidly decaying 
early light curves, e.g.. lNousek et al.ll2006l : IZhang et al. 1120061 ). In some cases, the afterglows 
are fainter at early times (t ~ 10^ — lO^s) than previously suspected. The explanation for 
why we see lines possibly correlated with flare emission may simply be that both depend on 
similar circumstances — non overpowering external shock emission — for their detection. 
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Fig. 1 6. — Equivalent width EW versus time in the source frame, adapted from iButler. et al. 



(J2005bl ). The four bursts studied here are marked as solid squares. The dotted curves show how 
the point for GRBs 060202, 050822, and 050714B move as the redshift is varied (small solid squares 
mark z = 1 for each burst). A number of points are measured for GRB 060218 in the ~16,000 
cts spectra (connected by a line), with one possible larger EW value from the ~ 500 cts spectra 
(Section U]). Previous detections and upper limits a re marked with open circles and la beled. The 
solid curves repr esent photoionization models from iBallantyne &: Ramirez-Ruizj (|200ll ) , explained 
in more detail in lButler. et al.l (|2005bl ). 
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Fig. 17. — Equivalent width EW versus time in the source frame for all (marginal) triggers at > 2a 
significance for 28 SwiftXI{T afterglows with measured redshifts, excluding GRB 060218. Redshifts 
are taken from the GCN. We restrict to only the spectra observed with exposures dt < tmean/2- 
Ninety percent of the EW values below 0.1, 0.01, or 0.01 days have EW < 1 keV. There are few 
spectra (3% of the sample) at times > 0.1 day in the source frame. 



Both the trigger times and apparent transient nature of the lines limit the distance from 
the progenitor: 

Ciobs 1 ^ 1 ril3 



R< 



1 + zl- cos(6') 



<;, 10^ 



cm, 



(2) 



wher e 6 is the angle between the line emitting material and the line of sight (e.g.. lVietri et al. 
20011 ). This implicates models from the so calle d "nearby reprocessor" class, whereby the 
lines are produced by a long-lasting central engine (JRees fc Meszarod 20001) or with the help o f 



magnetic field energy stored in the plasma bubble discussed above (JMeszaros fc Rees 



The small distaiices fro m the progenitor favor one-step "hypernova" explosion s fe.g..lWooslev 



1993 



2001). 



Paczynskil Il998l ) and disfavor a SN occurring prior to the GRB (e.g, IVietri fc Stella 
19981 ). For the plasma bubble, a possibly large magnetic field (~ 10^ Gauss) can be entrained 



and carried with the flow, both maintaini ng clumps of matter i n the flow and allowing 
for synchrotron irradiation of the clumps ( Meszaros fc Reesll200ll ). Or the lines could be 
generated by reflecti on from cold matter in the funnel wall carved out by the GRB jet from 
the progenitor (e.g.. iKallman. Meszaros. fc Reesll2003l ). 



The chemical species we infer p oint to ionization parameters ^ = L^/nR^ ^ 10^ 
(JLazzati. Ramirez- Ruiz, fc Reed |2002| ) . If the continuum luminosities Lx we measure are 



directly responsible for the lines, this implies a density n ~ 10 cm (GRB 060218) - 10 



^2l 
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cm~^ (GRB 050822). B ecause the largest dens i ties expected in the nea rby reprocessor sce- 
nario are ^ 10^^ cm"^ (JRees &: Meszarod l2000l : iMeszaros fc Reesll200ll ). a small fraction of 



the observed continuum (~ 0.01 — 10%) must be powering the lines. 



5.4. Systematic Uncertainties in the XRT Response 



Calibration efforts for the XRT spectral response are ongoing^. The quality of the 
current response matrices at low energy could have an important impact on the explanation 
for the soft excesses, particularly on the possibility of X-ray lines. Here we address a number 
of concerns raised by the XRT team (D. Burrows, private communication). This can be done 
quantitatively now due to the release of the v008 response files. The v008 files allow for the 
fitting of both the PC and WT mode data in the full 0.3-10 keV band, with the requirement 
of a small (3%) systematic error component. 

We have rerun the analysis using the v008 response files and including the systematic 
errors. For the 0.3-10 keV band, we find consistent distributions for the continuum fit 
parameters and continuum fit goodnesses (Figure [2]) and line trigger significances (Figure [3]) 
for the PC mode data. The distributions for continuum fit parameters for the WT mode 
data are also consistent. There is, however, a dramatic increase in the number of line triggers 
for the WT mode data for the well exposed spectra of 5 bright afterglows (GRBs 050502B, 
060124, 060202, 060210, 060211A). The line triggers fall near 0.5 keV and are clearly due to 
a dip near the 0-K edge at 0.532 keV, which can also been see in Figures [T2] and [HI for the 
spectra for GRBs 050714B and 060218. The XRT team argues that the dip — also seen in 
calibration observations of the Crab, 3C273, and H1426-I-428 — is due to a gain offset ^ 80 
eV, which occurs occasionally for WT mode data due to a problem with the bias subtraction. 
(It may also occur for the PC mode data, possibly due to illumination by the bright Earth.) 
Because we observe the magnitude of the dip to be proportional to the absorbing column 
Nh (hence to the depth of the 0-K edge), we agree with this explanation. This observation 
rules out the more pernicious possibility of a significant problem with the response files below 
0.6 keV. Also, we note that the spectra for the 5 WT mode bursts showing the dip can be 
adequately fit by absorbed power-laws if we also fit for the gain offset. 

One possible additional concern with 3 of the 4 bursts with sigificant line triggers is 
that the triggers appear to occur near the same observed energies (i.e., at ~100 eV intervals 
near 0.8, 0.9, 1.0 keV, etc), possibly indicating an instrumental origin. From simulations, 
we note that line-like residuals due to the gain offset are unimportant ^10% effects in this 



"'http://swift.gsfc.nasa.gov/docs/swift/analysis/xrt_digest.html 
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energy range. (There can be narrow residuals near the Si-K edge at 1.839 keV at the ~ 
20% level, however). Also, fitting for the gain, we find that the inferred trigger significances 
do not change. We, therefore, rule out the possibility that the instrumental effect of the 
gain offset produces the line triggers. Reinforcing this conclusion, a large number spectra 
with comparable source intensities to those with possible line detection are well fit by simple 
power-laws (Figures [2JA., [3]). However, we cannot entirely rule out the possibility that the low 
spectral resolution of the detector could allow for acceptable fits of an additive continuum 
component by closely spaced lines. This possibility is strengthened by the degradation of the 
of spectral resolution (from 0.08 keV at launch to ~ 0.1 keV in late 2006) due to radiation 
damage, an effect which is not treated in the v007 or v008 response files. 



6. Conclusions 

We have conducted a thorough and blind search for emission lines through nearly 0.4 
years of X-ray afterglow data accumulated by the Swift XKT. The majority of spectra (90%) 
are well fit by absorbed power-law models, but significant outliers to the population at the 5- 
10% level exist and have anomalously soft, possibly thermal spectra. Four bursts are singled 
out as possibly exhibiting 2-5 emission lines: GRBs 060218, 060202, 050822, and 050714B. 
Removed from the sample, line triggers in the spectra of the other bursts are consistent with 
Poisson fluctuations in the continua. 

The most significant soft component detections in the full data set of ~ 2000 spectra 
correspond to GRB 060218/SN 2006aj, with line triggers ranging from 4a to ~ 20cr. A 
thermal plasma model fit to the data indicates that the emission may be primarily due to 
L-shell transitions of Fe at ~ solar abundance. The possible line emission occurs near t ~ 1 
ksec, with a similar observed duration, indicating emitting material at i? ~ 10^^ cm. We 
suggest that the emission is due to the mildly relativistic cocoon of matter surrounding the 
GRB jet as it penetrates and exits the surface of the progenitor star. We associate the (> 4cr 
significant) line emission from 3 other events lacking redshifts with K-shell transitions in 
light metals. The lines in these bursts point to emission at similar distances, possibly at 
similar densities ~ 10^^ cm~^, and possibly subject to similar fluxes of ionizing radiation. 

As an alternate possibility — difficult to distinguish with the broad XRT spectral res- 
olution (~ 80 eV FWHM at 1 keV) — we successfully model the spectra using blackbody 
continuum components in addition to power-laws. With radii (again ~ 10^^ — 10^^ cm) and 
temperatures (~ 0.1 — 0.5 keV) implied by the fits, the possibility exists that the emission 
is continuous rather than discrete and possibly again due break out of the GRB shock or 
plasma cocoon from the progenitor star. We find that the energetics of the GRB and its 



-30- 



shock, inferred from the 4 events where breakout emission may be present, point toward 
the cocoon as the hkely source of the soft component. Bursts faint in 7— rays with spectra 
dominated by a soft X-ray flux possibly due to the shock breakout may outnumber classi- 
cal GRBs 100-1. Typically, the breakout flux would be faint or dominated by the external 
shock afterglow emission. The degeneracy between continuum and discrete emission compo- 
nents could possibly be lifted if redshifts are determined for the events discussed above or 
if redshifts are measured for bursts detected by Swift in the coming years with soft X-ray 
anomalies. 

N. Butler gratefully acknowledges support from a Townes Fellowship at U. C Berkeley 
Space Sciences Laboratory and partial support from J. Bloom and A. Filippenko. Special 
thanks to J. Bloom and the U. C Berkeley GRB team for comments on the manuscript 
and several useful conversations. Additional thanks to the Swift team for impressively rapid 
public release and analysis of the XRT data. 
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Table 1: Time-resolved Spectroscopy in ~ 16, 000 cts Spectra for GRB 060218 



Time (kscc) 



kT (kcV) 



1F7^ 



~7~r^ 



WT Mode 



1 


0.161-0.475 


211.00 ± 1.72 


0.94 ± 0.01 


1718.48/535 


372.40 ± 3.61 


1.54 ± 0.02 


3.97 ± 0.10 


595.99/535 






124.30 ± 53.90 


0.18 


+ 


324.80 ± 7.37 


1.50 ± 0.04 


5.73 ± 0.54 


537.20/530 


2 


0.475-0.688 


40.28 ± 2.71 


0.98 ± 0.01 


1662.32/549 


571.20 ± 5.47 


1.46 ± 0.02 


3.72 ± 0.10 


573.38/549 






147.18 ± 70.86 


0.17 


-1- 


496.20 ± 10.19 


1.44 ± 0.04 


5.31 ± 0.54 


535.00/544 


3 


0.688-0.872 


45.30 ± 3.06 


0.95 ± 0.01 


1724.45/529 


661.90 ± 6.35 


1.57 ± 0.02 


4.11 ± 0.10 


629.89/529 






340.95 ± 152.11 


0.16 


+ 


591.50 ± 13.57 


1.55 ± 0.04 


6.21 ± 0.56 


537.88/524 


4 


0.872-1.046 


48.00 ± 3.22 


0.92 ± 0.01 


1499.49/522 


689.20 ± 6.76 


1.62 ± 0.02 


3.97 ± 0.10 


640.10/522 






761.00 ± 348.78 


0.13 


+ 


659.20 ± 16.94 


1.71 ± 0.04 


6.85 ± 0.59 


603.46/517 


5 


1.046-1.223 


44.07 ± 2.96 


0.87 ± 0.01 


1444.86/500 


670.20 ± 6.87 


1.80 ± 0.02 


4.49 ± 0.10 


621.45/500 






2067.00 ± 851.39 


0.12 


+ 


691.80 ± 20.88 


1.94 ± 0.04 


8.26 ± 0.56 


539.78/495 


6 


1.223-1.411 


39.48 ± 2.65 


0.82 ± 0.01 


1486.17/493 


622.80 ± 6.66 


1.91 ± 0.02 


4.46 ± 0.10 


602.29/493 






1992.59 ± 849.83 


0.12 


-f 


652.90 ± 21.46 


2.03 ± 0.04 


8.04 ± 0.57 


529.29/488 


7 


1.411-1.617 


33.38 ± 2.23 


0.77 ± 0.01 


1419.64/470 


545.10 ± 6.16 


2.03 ± 0.02 


4.47 ± 0.09 


584.09/470 






6199.44 ± 2325.48 


0.10 


+ 


640.30 ± 22.94 


2.25 ± 0.04 


9.24 ± 0.50 


458.65/465 


8 


1.617-1.851 


27.12 ± 1.83 


0.71 ± 0.00 


1413.91/440 


478.90 ± 6.18 


2.23 ± 0.02 


4.64 ± 0.09 


566.43/440 






4650.94 ± 1842.21 


0.11 


+ 


568.40 ± 23.49 


2.42 ± 0.04 


8.83 ± 0.52 


446.16/435 


9 


1.851-2.116 


22.03 ± 1.48 


0.64 ± 0.00 


1563.68/422 


406.90 ± 5.56 


2.34 ± 0.02 


4.33 ± 0.09 


578.89/422 






4248.35 ± 1666.76 


0.11 


+ 


483.30 ± 21.47 


2.50 ± 0.04 


8.59 ± 0.51 


442.47/417 


10 


2.116-2.401 


18.26 ± 1.30 


0.57 ± 0.00 


1845.22/389 


380.70 ± 6.16 


2.57 ± 0.03 


4.45 ± 0.09 


640.99/389 






3146.32 ± 1199.94 


0.11 


+ 


417.70 ± 20.08 


2.61 ± 0.05 


8.10 ± 0.48 


400.39/383 


11 


2.401-2.753 


14.43 ± 0.98 


0.53 ± 0.00 


1791.20/371 


313.00 ± 5.57 


2.69 ± 0.03 


4.34 ± 0.08 


630.66/371 






2706.92 ± 986.90 


0.11 


+ 


337.90 ± 16.90 


2.71 ± 0.05 


7.88 ± 0.45 


344.64/365 



PC Mode 



12 5.950-2.872 Msec 



0.009 ± 0.001 
0.43 ± 0.16 



0.36 ± 0.02 
0.10 ± 0.01 



215.47/93 



0.031 ± 0.003 
0.03 ± 0.01 



3.41 ± 0.13 
3.17 ± 0.24 



4.39 ± 0.27 
7.33 ± 1.36 



144.03/93 
87.94/91 



Notes: Column density Nu measured in units of lO'^ 
(|Dickev fc LockmanlEoOol) . Nh fixed at lower limit 3.7 x 10^ 



iVff.Galactic = 1.11 X 10^ 

1^^ for pure blackbody (BB) fits. Flux 
measured in units of 10^^^ erg cm^^ s^^, 0.5-10 keV for the power-law (PL) and bolometric for the black- 



body. Temperature measured in the source frame, z = 0.033 ( Mirabal et al.ll2006r ). 
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Table 2: Time-resolved Line Search in ~ 16, 000 cts Spectra for GRB 060218 



# Time N Signif. Line Energy , Equivalent Width , Flux 

[s] lines (keV,eV,10~^'^ erg cni~^ s~^) 

WT Mode 



1 


0.161-0.475 





<5.0ct 


2 


0.475-0.688 





< 5.0cr 


3 


0.688-0.872 





< 5.0(7 


4 


0.872-1.046 





<5.0cr 


5 


1.046-1.224 


5 


7.0(7 


6 


1.224-1.411 


5 


5.9(7 


7 


1.411-1.618 


5 


8.2(7 


8 


1.618-1.851 


5 


8.0(7 


9 


1.851-2.116 


5 


9.3(7 


10 


2.116-2.401 


5 


13.3(7 



(0.93,28,5.5) (3.19,45,3.3) (3.65,40,2.6) (1.02,20,3.6) (2.76,30,2.5) 
(0.95,26,5.2) (1.14,16,2.8) (3.24,28,1.9) (0.87,23,4.8) (1.03,16,3.0) 
(0.89,42,9.1) (0.80,40,9.7) (1.06,22,3.9) (4.04,47,2.1) (0.97,18,3.5) 

(1.00,37,7.3) (0.79,51,13.4) (0.90,35,7.9) (1.12,23,3.9) (0.71,47,14.1) 
(0.96,34,6.0) (0.88,32,6.4) (0.79,33,7.5) (1.06,21,3.3) (2.55,31,1.5) 

(0.92,55,10.7) (1.01,40,6.8) (0.80,39,9.3) (1.11,26,3.8) (0.75,33,8.8) 



11 


2.401-2.7535 


5 


14.1(7 


(0.90,49,8.4) (0.99,40,5.8) (0.78,42,9.1) (1.12,23,2.8) (0.85,24,4.6) 


PC Mode 


12 


5.950-2.872 Msec 


4 


7.1(7 


(0.80,109,0.003) (0.89,78,0.001) (0.69,96,0.003) (1.07,43,0.001) 
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Table 3: Time-Integrated Spectroscopy for GRB 060218 



Model 



Nh 



Norm. 



r or kT (keV) 



xV^ 



WT Mode Data t = 160s-2.8ksec 



power-law 
blackbody 
power-law 

power-law 



blackbody 
MEKAL 



3.6 ± 0.1 
0.4 ± 0.1 
7.2 ± 0.2 

4.0 ± 0.2 



(3.43±0.01) xlO-9 
(1.84±0.01 xlO^^ 

(3.51±0.03 xlO-9 
(7.6± 0.9) xlO-9 

(3.40± 0.01) xlO'^ 
(8.0± 0.5) xlO-2 



1.84 ± 0.01 
0.769 ± 0.002 

1.95 ± 0.01 
0.123 ± 0.002 

1.82 ± 0.01 

0.82 ± 0.02 



1740.56/768 
14434.53/768 

1111.33/766 

1323.00/766 



PC Mode Data t = 5.9ksec-3.1Msec 



power-law 
blackbody 
power-law 

power-law 



blackbody 
MEKAL 



4.4 ± 0.3 
0.5 ± 0.2 
6.9 ± 1.3 

5.6 ± 0.7 



(3.1±0.3) xlO-13 
(9.0±0.4) xlO~" 
(2.5±0.6) xlO-13 
(2.9±0.7) xlO-12 
(2.5±0.32) xlO-" 
(2.8±0.5) xlO~4 



3.4 ± 0.1 
0.36 ± 0.02 

3.1 ± 0.2 
0.10 ± 0.01 

3.1 ± 0.1 
0.23 ± 0.03 



143.99/93 
215.24/93 

87.63/91 

66.01/91 



Notes: Column density Nh measured in units of 10^^ cm^^. -/Vff.caiactic = l-H x 10^^ cm^^ 
( Dickey &: LockmanI 119901) . Power-law model normalization (Norm.) in units of erg cm^^ s^-'^ [0.5-10 
keV]. Blackbody temperature measured in the source frame. Blackbody normalization in units of erg cm^^ 



MEKAL ( Mewe. Gronenschild. fc van den Oordlll985[) model normalization in units of and equal to 



i^TDjir+iy 



■ / neUHdV , for angular diameter distance Da, electron density Ue, and proton density uh- 
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Table 4: Time-resolved Spectroscopy for GRB 050822 



# 


Time (s) 


BB Flux 


kT (keV) 


Nh 


Tl 


PL Flux 


r 


Nh 


x'/^ 


WT Mode 


1 


111.0-121.0 


26.1 ± 1.5 


0.66 ± 0.03 


0.1 ± 0.1 


105.93/59 


44.6 ± 3.4 


1.8 ± 0.2 


1.9 ± 0.6 


48.50/59 


2 


121.0-131.0 


28.0 ± 1.6 


0.63 ± 0.03 


0.1 ± 0.1 


107.13/66 


51.1 ± 4.0 


2.0 ± 0.2 


2.7 ± 0.6 


53.64/66 


3 


131.0-141.0 


26.1 ± 1.4 


0.57 ± 0.03 


0.1 ± 0.1 


134.80/65 


50.0 ± 4.2 


2.1 ± 0.2 


2.5 ± 0.6 


64.32/65 


4 


141.0-151.0 


20.5 ± 1.3 


0.54 ± 0.03 


0.1 ± 0.1 


103.56/54 


38.6 ± 3.2 


2.0 ± 0.2 


1.9 ± 0.6 


32.66/54 


5 


151.0-160.9 


16.5 ± 1.2 


0.47 ± 0.03 


0.1 ± 0.1 


85.38/46 


33.6 ± 3.9 


2.4 ± 0.3 


2.4 ± 0.7 


33.79/46 


6 


160.9-170.9 


10.1 ± 0.8 


0.40 ± 0.04 


0.1 ± 0.1 


80.37/31 


18.2 ± 2.6 


2.2 ± 0.4 


1.2 ± 0.9 


44.87/31 


7 


170.9-200.9 


3.8 ± 0.3 


0.35 ± 0.03 


0.1 ± 0.1 


91.20/36 


6.4 ± 0.8 


2.3 ± 0.3 


0.7 ± 0.7 


33.05/36 


8 


200.9-240.8 


3.6 ± 0.3 


0.37 ± 0.02 


0.1 ± 0.1 


80.67/44 


7.3 ± 1.2 


2.7 ± 0.3 


2.2 ± 0.7 


39.59/44 


9 


240.8-280.8 


2.7 ± 0.2 


0.35 ± 0.03 


0.1 ± 0.1 


47.69/32 


5.5 ± 1.3 


3.1 ± 0.4 


2.2 ± 0.8 


24.38/32 


10 


280.8-429.6 


2.1 ± 0.2 


0.34 ± 0.02 


0.1 ± 0.1 


76.08/52 


4.6 ± 1.0 


3.1 ± 0.3 


2.4 ± 0.7 


51.37/52 


11 


429.6-449.6 


10.3 ± 2.3 


0.26 ± 0.02 


0.7 ± 0.6 


33.37/43 


81.1 ± 44.3 


5.2 ± 0.5 


5.6 ± 1.0 


30.33/43 


12 


449.6-469.5 


10.0 ± 2.0 


0.25 ± 0.02 


0.4 ± 0.5 


35.46/46 


99.6 ± 54.3 


5.6 ± 0.5 


5.7 ± 1.0 


41.03/46 


13 


469.5-489.5 


10.2 ± 4.1 


0.20 ± 0.02 


0.8 ± 0.8 


35.17/37 


169.7 ± 72.5 


6.7 ± 0.8 


6.6 ± 1.4 


35.03/37 


14 


489.5-509.4 


6.4 ± 0.7 


0.21 ± 0.02 


0.1 ± 0.4 


51.06/33 


53.3 ± 11.1 


5.9 ± 0.9 


4.8 ± 1.6 


54.71/33 


15 


509.4-549.3 


4.1 ± 0.4 


0.17 ± 0.02 


0.3 ± 0.6 


31.63/31 


23.6 ± 6.6 


6.1 ± 1.0 


4.2 ± 1.5 


37.11/31 


PC Mode 


16 


354.4-12823 


0.45 ± 0.03 


0.27 ± 0.02 


0.1 ± 0.1 


145.96/38 


0.81 ± 0.09 


2.1 ± 0.2 


0.6 ± 0.2 


52.01/38 


17 


12823-34706 


0.19 ± 0.02 


0.39 ± 0.05 


0.1 ± 0.1 


107.85/34 


0.40 ± 0.04 


2.1 ± 0.3 


1.3 ± 0.6 


45.42/34 


18 


34706-88110 


0.04 ± 0.01 


0.43 ± 0.04 


0.1 ± 0.1 


72.60/32 


0.08 ± 0.01 


2.2 ± 0.3 


1.6 ± 0.6 


31.15/32 


PC + WT 


19 


549.3-1000.0 


1.2 ± 0.3 


0.17 ± 0.02 


0.3 ± 0.1 


48.84/19 


12.6 ± 1.3 


4.0 ± 1.1 


1.2 ± 1.4 


35.12/19 



Notes: Column density Nh measured in units of 10^ 



iV, 



//, Galactic 



^ 2.25 X 10^ 



( Dickey &: LockmanlllQQOl ). Power-law (PL) and blackbody (BB) fluxes measured in 10 ^° erg 
0.5-10 keV band for PL model. Blackbody temperature measured in the observer frame. 
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Table 5: Time-resolved Line Search for GRB 050822 



# 



Time 

[si 



N Signif. 
lines 



Line Energy , Equivalent Width , Flux 
(kcV,eV,10-" erg cm-^ s'^) 



WT Mode 



1 


111.0-121.0 





< 1.8(7 


2 


121.0-131.0 





< 1.7cr 


3 


131.0-141.0 


3 


2.1cr 


4 


141.0-151.0 





< 1.2cr 


5 


151.0-160.9 


2 


2. Oct 


6 


160.9-170.9 


4 


2.7cr 


7 


170.9-200.9 





< 1.9cr 


8 


200.9-240.8 





< 1.5o- 


9 


240.8-280.8 





< 0.6cr 


10 


280.8-429.6 





< 1.20- 


11 


429.6-449.6 





< 0.6cr 


12 


449.6-469.5 





< 1.3cr 


13 


469.5-489.5 





< 1.6cr 


14 


489.5-509.4 


5 


4.4cr 


15 


509.4-549.3 





< 1.9cr 



(1.94,93,7.7) (2.91,115,5.9) (3.66,192,7.6) 

(1.01,49,5.9) (1.11,51,5.5) 
(0.66,90,7.6) (0.88,67,4.2) (1.13,111,5.5) (1.29,95,4.1) 



(0.81,82,3.6) (0.91,142,4.5) (1.04,194,4.3) (1.23,221,3.1) (1.49,265,2.2) 



PC Mode 


16 

17 
18 


354.4-12823 
12823-34706 
34706-88110 


3 





2.5cr 

< 1.9cr 

<2.0f7 


(0.41,75,0.6) (0.63,106,0.5) (3.48,457,0.3) 


PC + WT 


19 


549.3-1000.0 


5 


3. Oct 


(0.65,131,1.1) (0.75,64,0.4) (0.82,65,0.3) (0.90,65,0.3) (1.09,50,0.2) 



Table 6: Time-resolved Spectroscopy for GRB 050714B 



Time Reg. (s) 



feT (kcV) 



Nh 



xV" 



Nh 



X^/" 



WT Mode 



157.4-189.7 



7.1 ± 1.! 



0.23 ± 0.02 2.66 ± 0.59 33.49/31 85.49 ± 40.36 5.90 ± 0.44 0.81 ± 0.09 29.39/31 



PC Mode 



2 221. 4-58. 3ks 0.006 ± 0.001 0.30 ± 0.03 0.30 ± 0.23 97.19/35 0.015 ± 0.003 3.24 ± 0.28 0.31 ± 0.05 76.78/35 

0.033 ± 0.007 0.13 ± 0.03 3.8 ± 1.4 -I- 0.008 ± 0.002 2.2 ± 0.5 tied 48.64/33 

3 275.0-525.0 0.74 ± 0.37 0.22 ± 0.03 0.83 ± 0.69 35.67/19 4.28 ± 2.69 5.11 ± 0.63 0.53 ± 0.10 39.92/19 
7.8 ± 0.4 0.11 ± 0.03 4.5 ± 2.3 + 0.7 ± 0.1 3.1 ± 1.1 tied 21.40/17 



Notes: Column density Nh measured in units of 10^^ cm ^. iV/f, Galactic = 5.31 x 10"' 
( Dickey fc Lockmanlll990l ). Flux measured in units of 10^^° erg cm^'^ s^^, 0.5-10 keV for the power-law 



(PL) and bolometric for the blackbody (BB). Blackbody temperature measured in the observer frame. 
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Table 7: Time-resolved Line Search for GRB 050714B 



# 



Time 

[si 



N 
lines 



Signif. 



Line Energy , Equivalent Width , Flux 
(keV,eV,10-i2 erg cm-^ s"^) 



WT Mode 



157.4-189.7 



<1.8cr 



PC Mode 



2 221.4-58.3 ks 4 4.2cr (0.91,147,0.09) (0.76,141,0.12) (0.56,238,0.35) (1.12,64,0.03) 

3 275.0-525.0 4 3.4cr (0.74,170,57.5) (0.59,340,261) (0.93,108,16.0) (0.83,68,15.4) 



Table 8: Spectroscopy of Bright (> 10 cts/s peak) XRT Flares 



Burst 


Time Reg. (s) 


Flucncc 
(10^^ erg cm^^) 


r 


Nh 
(1021 cm-2) 


2 / 


050502b 


400-1200 


6.61 ± 0.09 


2.43 ± 0.03 


1.76 ± 0.07 


373.28/357 


050607 


250-600 


0.38 ± 0.03 


2.29 ± 0.13 


2.11 ± 0.36 


30.78/34 


050712 


150-300 


0.55 ± 0.02 


2.06 ±0.10 


2.35 ± 0.32 


91.62/73 


050713A 


95-150 


3.07 ± 0.10 


2.31 ± 0.06 


5.55 ± 0.25 


216.45/207 


050714B 


275-525 


1.07 ± 0.67 


5.11 ±0.63 


5.25 ± 1.02 


39.92/19 


050730 


130-300 


1.43 ± 0.04 


1.63 ± 0.05 


1.69 ± 0.21 


178.00/171 


050730 


300-600 


1.90 ± 0.04 


1.72 ± 0.05 


1.02 ± 0.14 


215.10/227 


050730 


600-800 


0.84 ± 0.03 


1.84 ± 0.07 


0.43 ± 0.18 


127.01/134 


050820 


215-252 


2.40 ± 0.06 


0.99 ± 0.05 


1.87 ± 0.23 


202.11/207 


050822 


410-650 


4.03 ± 0.50 


4.95 ±0.13 


4.07 ± 0.24 


141.36/127 


050904 


350-600 


1.37 ± 0.03 


1.73 ± 0.05 


1.49 ± 0.15 


221.51/227 


051117A 


1250-1725 


6.81 ± 0.09 


2.31 ± 0.03 


1.94 ± 0.08 


368.05/364 


051117A 


800-1250 


4.77 ± 0.08 


2.32 ± 0.03 


1.76 ± 0.09 


333.41/307 


060111A 


200-500 


5.59 ± 0.09 


2.29 ± 0.03 


2.70 ± 0.10 


355.60/335 


060124 


300-650 


16.48 ± 0.13 


1.39 ± 0.02 


2.36 ± 0.07 


611.60/611 


060124 


650-900 


12.91 ± 0.13 


1.92 ± 0.02 


2.26 ± 0.08 


486.54/453 


060204B 


100-270 


1.95 ± 0.05 


2.00 ± 0.05 


2.78 ± 0.19 


168.65/201 


060204B 


270-450 


0.25 ± 0.02 


2.68 ±0.13 


2.88 ± 0.34 


78.40/70 


060210 


100-165 


1.32 ± 0.03 


1.82 ± 0.06 


2.04 ± 0.20 


149.31/157 


060210 


165-300 


2.70 ± 0.05 


1.92 ± 0.04 


1.98 ± 0.14 


210.08/260 


060210 


350-450 


1.18 ± 0.07 


2.93 ± 0.08 


2.86 ± 0.21 


127.92/141 


060312 


100-200 


0.14 ± 0.01 


2.30 ± 0.09 


2.99 ± 0.29 


72.11/105 
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Table 9: Line Detections in the Bright (> 10 cts/s peak) XRT Flares 



Burst 



Time 

[si 



N 
lines 



Signif. 



Line Energy , Equivalent Width , Flux 
(keV,eV,10~^^ erg cm"^ s~^) 



050502b 


400-1200 


4 


3.6(7 


050607 


250-600 


5 


2.1(7 


050712 


150-300 





<1.9(7 


050713A 


95-150 


3 


2.9(7 


050714B 


275-525 


3 


3.4(7 


050730 


130-300 


1 


2.2(7 


050730 


300-600 


5 


2.9(7 


050730 


600-800 





<1.9(7 


050820 


215-252 


2 


2.3(7 


050822 


410-650 


3 


3.9(7 


050904 


350-600 





<1.7a 


051117A 


1250-1725 


3 


4.4(7 


051117A 


800-1250 


4 


4.3(7 


060111A 


200-500 


4 


2.3(7 


060124 


300-650 


3 


3.2(7 


060124 


650-900 


5 


2.6(7 


060204B 


100-270 


4 


2.8(7 


060204B 


270-450 





<1.8(7 


060210 


100-165 





<1.6(7 


060210 


165-300 





<1.3(7 


060210 


350-450 


1 


2.1(7 


060312 


100-200 


1 


2.5(7 



(1.44,15,0.4) (1.69,18,0.4) (4.48,57,0.3) (2.14,17,0.2) 
(2.11,152,0.2) (0.79,62,0.4) (0.61,118,1.1) (3.29,409,0.3) (2.57,216,0.3) 

(1.82.40.4.1) (4.48,97,3.0) (3.72,65,2.6) 
(0.74,162,6.9) (0.59,301,30.7) (0.93,98,1.7) 

(1.97,53,0.7) 
(1.12,31,0.4) (1.41,29,0.3) (1.01,27,0.4) (1.87,47,0.5) (0.77,18,0.3) 

(4.27,72,4.9) (1.96,41,2.8) 
(0.91,22,2.8) (0.61,44,30.9) (2.06,76,0.3) 

(1.88,34,0.9) (2.83,26,0.4) (0.60,25,3.0) 
(1.75,30,0.6) (1.37,20,0.5) (0.92,16,0.7) (1.05,12,0.4) 
(1.80,20,0.7) (1.04,12,0.9) (0.90,11,0.9) (0.80,10,1.0) 

(1.89.18.1.2) (4.57,29,1.4) (3.47,18,0.9) 
(3.20,28,1.5) (4.51,40,1.6) (4.89,42,1.6) (1.12,8,1.2) (2.73,18,1.1) 

(2.64,61,1.7) (2.08,44,1.5) (2.83,56,1.4) (4.32,83,1.4) 



(0.62,50,8.7) 
(0.60,117,1.7) 



Table 10: Afterglow Energetics Parameters 



Burst 



-^7,iso 

[10^1 erg] 



Lx,t=10hrs SbB 

[10^'' erg s-i] [10^1 erg] 



060218 0.033 

060202 0.783 
050822 1.2 

050714B 2.66 



0.062 ±0.003 

4.4 ±0.4 

15 ±1 

13 ±2 



0.013 ±0.003 
43 ±11 
410 ± 80 
440 ± 90 



StBB[s] 

[b] 



0.23 ±0.07 -300 

1.3 ±0.4 j^lOOO 

0.43 ±0.18 -100 

2.8 ±1.0 -200 



Notes: Redshifts z for GRBs 060202, 050822, and 050714B are inferred from the X-ray spectroscopy (Section 
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